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Abstract 

Supersymmetric Yang Mills theory is directly accessible to lattice simulations using 
current methodology, and can provide a non-trivial check of recent exact results in 
SQCD. In order to tune the lattice simulation to the supersymmetric point it is nec- 
cessary to understand the behavior of the theory with a small supersymmetry breaking 
gaugino mass. We introduce a soft breaking gaugino mass in a controlled fashion using 
a spurion analysis. We compute the gluino condensate, vacuum energy and bound-state 
masses as a function of the gaugino mass, providing more readily accessible predictions 
which still test the supersymmetric results. Finally we discuss diagnostics for obtaining 
the bare lattice parameters that correspond to the supersymmetric continuum limit. 
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1 Introduction 



The recently proposed exact solutions ^ of N=l supersymmetric QCD (SQCD) satisfy 
striking tests of self-consistency and provide an extremely plausible picture of the rich low- 
energy dynamics of these models. Nevertheless, one may feel a little discomfort at the 
absence of direct non-perturbative tests of the results. An obvious possibility is that these 
models could be simulated directly on the lattice. Some initial work in these directions 
has already been performed in [Q. Unfortunately, as is well known, lattice regularization 
violates supersymmetry 0, and a special fine-tuning is required to recover the SUSY limit 
(this is analogous to the case of chiral symmetry in lattice QCD). Away from the SUSY 
point, the continuum limit of the lattice theory is described by a model with explicit SUSY 
violating interactions. In some cases, these violations may correspond only to soft breakings 
0, although this is not guaranteed in general. 

Softly broken SUSY models can be studied using spurion techniques, and "exact" re- 
sults are possible || (additional investigations of softly broken SQCD can be found in 
0). In this paper, we propose some tests of supersymmetric Yang-Mills theory (SYM - or 
SQCD with zero flavors of quark) which can be carried out using lattice techniques. SYM 
is a simple theory with only one parameter, the gauge coupling. The only low-dimension 
(renormalizable) SUSY violation allowed by gauge invariance is a gaugino mass, which is a 
soft violation. Therefore, the continuum limit of the lattice regularized version of SYM is 
simply SYM with a massive gaugino. The SUSY limit can be reached by fine-tuning the 
lattice parameter corresponding to a bare gaugino mass. In order to understand this limit 
as well as possible, we study continuum SYM with explicit gaugino mass, and derive some 
relations describing the approach to the SUSY limit. The spurion techniques used assume 
that there is not a phase transition to some totally new phase the moment that supersym- 
metry is broken, a fact that may be tested on the lattice (if such a transition did occur then 
testing the supersymmetric predictions would become very hard since that phase would con- 
stitute only a single point in parameter space and there would be no understanding of the 
approach to that point). Assuming no such transition, several non-trivial predictions can 
be made regarding the vacuum energy and of the behavior of the gaugino condensate. A 
less rigorously derived description of the lightest bound states of SYM theory has also been 
proposed in the literature |jlO|, ^ from which predictions for the masses of the gluino-gluino 
and glue-gluino bound states and their splittings away from the supersymmetric point may 
be obtained. We also discuss aspects of tuning toward the SUSY limit. 

A lattice test of the pure glue theory would also provide a test of the tower of SQCD 
theories with fundamental matter flavors. The numerical coefficients in the effective super- 
potential of SQCD theory with Nf = Nc — 1 have been determined analytically |jl2| (this is 
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the first value of Nj wliere tlie gauge group can be broken completely and the theory studied 
in the perturbative regime). The pure glue theory predictions resulting from integrating out 
the quark flavors in that effective theory are closely related to these coefficients as we review 
below. 

The outline of the paper is as follows. In section 2 we review the holomorphic analysis of 
SYM. In sections 3 and 4 we introduce soft breakings into the analysis, deriving predictions 
for the lattice simulations. In section 5 we discuss diagnostics for tuning towards the SUSY 
point. In section 6 we summarize our results. Finally, in the appendix, we discuss the 
rescaling anomaly and the holomorphic vs. canonical field normalizations, which are relevant 
to the comparison between exact and lattice results. 



2 SUSY Yang-Mills 



The bare Lagrangian of SYM with SU{Nc) gauge group is 

^ ^ (1) 



^ = - 

9o 



This model possesses a discrete global Z2Nc symmetry, a residual non-anomalous subgroup of 
the anomalous chiral U{1). The theory is believed to generate a gaugino condensate breaking 
the Z2Nc symmetry to Z2 and also to exhibit a mass gap due to confinement. 

In supersymmetric notation the Lagrangian (|l|) can be written as 

£ = f (fe —ImroWWa , (2) 
J Stt 

where the gauge coupling is defined to be tq = ^ + 

In this notation, the strong coupling scale of SYM is defined through the relation 

A = e^'^'o/^oAt/y , (3) 

where 60 = SA^^c is the first coefficient of SYM /3-function and Auv is the UV cut off of the 
theory at which the coupling takes the value go- The relation of this definition of the strong 
scale to the MS scheme has been investigated in |]13[. Note that A^" is explicitly 2'7r-periodic 
in the 6'0-angle. For the purposes of this paper we set Oq = 0. 

To derive the low energy effective theory of SQCD we note that there are two anomalous 
symmetries of the theory, U{1)r and scale invariance. In fact their anomalies are related 
since their currents belong to the same supermultiplet. These symmetries can be restored 
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in an enlarged theory provided we allow the spurion couplings to transform: 

U{1)r : W{x, 6) e'^W{x, Oe'^) 

A ^ Ae^2C/3 



Scale Invariance : W{x, 6) e^^''^W{xe^, Oe^/^) 

A ^ Ae« 

We may now determine the general form of the partition function (assuming a mass gap) 
in the large volume limit as a function of r subject to these symmetries. The only possible 
terms are 



Z[t\ = exp iV 



-A^A\d- {pA^\F + h.c. 
a 



(4) 



The numerical coefficients a and f3 remain undetermined from the above symmetry ar- 
guments. f3 may be determined from the results for SQCD with massive quarks. The U{1)r 
symmetries of models with matter are sufficient to determine the form of the superpotential 

as 

det(m)A''«l^^^' (5) 



where A is defined similarly to (|^) but with the appropriate /5-function bo = 3Nc — Nf and 
c is a constant, undetermined by the symmetry. The scale transformation anomaly is not 
sufficient to uniquely determine the Kahler potential in these theories with two mass scales A 
and m. The coefficient c may be explicitly calculated in theories in which Nf > Nc — 1 since 
the squark vevs may be chosen to completely break the gauge symmetry in the perturbative 
regime. The superpotential term is then generated by calculable instanton effects. The 
calculation was performed in |]I2|. The mass terms may then be taken to infinity removing 



the quarks from the theory and leaving the pure glue theory in a controlled fashion. The 
resulting prediction for c in Minkowski space is {Nj — N^) and hence (3 is N^. 

These strong arguments lead to two predictions for the condensates of the SYM theory. 
The source J for the gaugino correlator AA occurs in the same position as the F-component 
of T and is hence known. There are two independent correlators 

(AA) = -327r2A3 
— 1024? 77 

(AAAA) = -^^IMW . (6) 

The IR theory has a gaugino condensate ~ A^, with phase 27rir/A^c and hence there are 
Nc degenerate vacua associated with the A^^cth roots of unity. Below, therefore, A^ is an Nc 
valued constant with phases n2ni/Nc where n runs from 0... — 1. 
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3 Soft Supersymmetry Breaking 



We may induce a bare gaugino mass through a non zero F-component of the bare gauge 
couphng t = tq + FrOe 

- ^/m[F.AA] (7) 

To make the gaugino mass canonical we take Fr = iSTcmx. In the IR theory r enters through 
the spurion, A^, the lowest component of which is the strong interaction scale A 



A^° = A^'o (1 - 16n^mx99) 



A occurs linearly in the superpotential of the theory. Thus there will be a correction to 
the potential of the form: 



AV 



a 



(9) 



Terms with super derivatives acting on the spurion field can also give rise to contributions to 
the potential but these are higher order in an expansion in mx/A. The shift in the potential 
energy of the A^^^ degenerate vacua of the SYM theory at linear order in mx is known and 
we may determine the vacuum structure 

"2?™ 



Al^ = -327r^|mAA^| cos 



(10) 



For small soft breakings, mx <^ A, where the linear term dominates, the degeneracy between 
the SYM vacua is broken favoring one vacuum dependent on the phase of the gaugino mass. 
The coefficient in the energy shift is a test of the exact superpotential in Eq.(H). 



We may also determine the leading shift in the gaugino condensate 

(AA) = 



327r2A3 + ^^m\\A\^ 



am 



(11) 



which depends on the unknown parameter a. Strictly speaking there are also divergent 
contributions to this quantity which are proportional to mx times the cut-off squared. We 
will have more to say about these divergences in section 5. 

Reinserting the bare Oq angle into the expression for the shift in vacuum energy we find 

27m 



AV = -327r^|mAA''| cos 



(12) 



As ^0 is changed first order phase transitions occur at 6q = (odd)7r where two of the SYM 
vacua interchange as the minimum of the softly broken theory. 



5 



4 The Lightest Bound States 



An alternative description of the low energy behaviour of SYM theory has been presented by 



Veneziano and Yankelowicz which attempts to describe the lightest bound states of the 
theory. The form of their effective action can be rigourously obtained from the discussion 
above. Since the source J for WW occurs in the same places as the coupling r we also know 
the source dependence of Z. If we wish we may Legendre transform Z[r, J] to obtain the 
effective potential for the classical field 

1 



S 



327r' 



■ Tr {W') 



(13) 



We find 

T[r,S] = ^[SS)"''\^ - N,[S-SHS/M'))\^+\i.c. . (14) 

So derived, this effective action contains no more information than Eq.(|^) simply being a 
classical potential whose minimum determines the vev of 5* and we find, by construction, 
Eq®. 

A stronger interpretation can also be given to the VY action. We can attempt to use 
it to describe the physical gluino-gluino and gluino-glue bound states, if we assert that 
they interpolate in the perturbative regime to the field WW. Under this assumption the 
symmetry properties of the bound states would be those of the S field, reproducing the VY 
action for those states, up to additional scale-invariant Kahler terms. (The latter only affect 
our predictions at higher order in the soft breaking.) To obtain the physical states one 
performs an appropriate rescaling of the S'-field 

3 



^1/3 



a 



(15) 



in the Lagrangian (|T^ to make the kinetic term canonical 



C 



D 
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:$3 _ ^3jj^ 



,a 



1/2 $ 



3 A' 



+ h.c. 



(16) 



In fact, as discussed in [0, this effective Lagrangian is not complete since it does not 
possess the full Z2Nc symmetry of the quantum theory. To restore that symmetry the extra 
term 



A/: 



27cim 



S-S 



(17) 



where m is an integer valued Lagrange multiplier must be added. For the n = vacuum 
with vanishing phase this extra term vanishes and the VY model above is recovered. We 
shall concentrate on that vacuum and real, positive mass perturbations which make that 
vacuum the true vacuum of the perturbed theory. 
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In addition to the particle states associated with the field (gluino-glunio or gluino- 
glue balls), one might expect to find states associated with the familar glueballs of QCD, 
with interpolating fields F"^ or FF, and their superpartners. A chiral supermultiplet can 
be constructed out of D^W^ = D^S which contains the appropriate fields. We have no 
reason to expect that the glueball states are parametrically heavier than the (f) field states, 
and so there is no a priori reason to integrate them out of the effective Lagrangian for 0. 
Furthermore, any glueballs are likely to be strongly coupled to and have a non-negligible 
effect on its dynamics. 

On the other hand, any couplings to glueballs are still constrained by scale invariance and 
anomalous f/(l)/j-invariance. Let us postulate the existence of some effective glueball chiral 
superfield 0'(^)) which is dimension one and has a canonical kinetic energy term. Since 0'(x) 
is constructed from D^W^ it must have R-charge zero. A superpotential term describing the 
interaction of 0' with S would be of the form: 



The constraint from anomalous R-invariance (under which A has charge 2/3), requires that 



or 6 = 0. Therefore there are no superpotential terms of the correct type. All interactions 
between the potential glueball excitations and S (including mass or wavefunction mixing) 
must appear as Kahler terms in the effective Lagrangian. There are already unknown Kahler 
terms that enter into the physical masses. These are higher derivative terms which, while 
suppressed by powers of A, contribute to the wave function renormalization, ZijP')^ which 
must be evaluated at = m? ~ A^ to obtain the physical masses. The Kahler corrections 
from the 0' may be subsumed into these unknown terms. 

Of course, it is possible that the glueball fields do not enter the effective Lagrangian 
as 0'(x), but rather as some other effective field, or even the auxilliary field of S (which 
contains G^j, and GG). In this case the description advocated here may not be correct. We 
will discuss this possibility (more specifically, the status of glueballs in SYM) in more detail 
below, but for now we simply adopt the VY Lagrangian and discuss its mass predictions. A 
lattice simulation will hopefully test these predictions and shed light on whether the action 
is indeed the correct description. 

The straightforward evaluation of bosonic (AA) and fermionic {gX) excitation masses 
around the minimum from Eq(|16|) gives 




(18) 



2 



2a-2/3(3a + 6-3) , 



(19) 



mx\ = mgx 



N^aA . 



(20) 



7 



It is important to stress again that these masses are not the physical masses of the bound 
states. Rather, they are zero-momentum quantities, which are related to the physical ones 
by wave function factors Zlp^ = fripf^yg). These wave function factors result from higher- 
derivative Kahler terms in C, and are unknown. 

A soft breaking gaugino mass may again be introduced through the F-component of the 



spurion (first investigated in fTTl)- The new scalar potential is 
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,a 



1/2 



327r2a3/2 



V = -^\r\ ln(^^) - — Remxr (21) 



27 



3 A' 

from which we can calculate the shifts in the masses of the bound states. The two scalar 
fields and the fermionic field are all split in mass. The eigenvalues of the mass matrices are 

167r^|mA| 



Mfermion = N^aA + 
Mscalar = A^c^A 



567r^|mA| 

21 



AOtt \m\ 

Mp.scalar = A^.aA + ^' (22) 

these results have been derived for real, positive mass which favor the supersymmetric vac- 
uum with vanishing phase. For even there is a supersymmetric vacuum characterized by 
(S) = — A'^ which is prefered by negative, real mass perturbations. It is easy to construct an 
effective lagrangian about that vacuum in a similar fashion to (0). The bound state masses 
are again given by (P^). 

The physical masses are again related to these quantities by unknown wave function 
renormalizations Z which arise from Kahler terms, 

^physical = Z M . 

Fortunately, we know that in the SUSY limit the wavefunction factors are common within 
a given multiplet. This degeneracy holds even after the vev of the field is shifted by the soft 
breakings since a shift in the vev alone (without SUSY breaking) leaves the physical masses 
degenerate within a multiplet. We also know that the relative change in these Kahler terms 
is of order 0{f^), and hence can be ignored at leading order in the soft breakings. Therefore, 
we may still obtain a prediction for the rate of change of the ratios of the physical masses, 

M^ra,) ^ Z^x)Mim,) - Z(0)M(0) 



Z(0)M(0) 

dM " 



1 1 dZ 
M ^ 'ZdM 



mx (23) 
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near the SUSY limit. The factor in brackets is common within a given multiplet. Since 
the quantity Z{mQ) is unknown, we can only predict the ratios of M at the SUSY point or 
equivalently the ratios of dM/ dm\ 

dMpJdmx : dUferm/dmx : dMjdmx = 5:6:7 (24) 

At this point it is worthwhile to return to the question of glueball masses and their 
interaction with the S bound states. We make an observation based on some results of 
West's on glueball masses in QCD Using QCD inequalties, West has shown that the 
mass of the lightest non-vacuum state coupling to the operator is less than the mass of 
the lightest non- vacuum state coupling to GG. In QCD, this implies that the lightest glueball 
is a scalar, not a pseudoscalar. West's results can be applied to SYM due to the positivity 
of the gluino determinant [|15[. In SYM the results are relevant both to the glueballs and 



the gluino-gluino bound states, as they have the same quantum numbers and can mix even 
at the perturbative level. Now, if ( P^ is correct, then very close to the SUSY point the 
shift in S bound state masses is such that the pseudo-scalar S bound state is lighter than 
the scalar S bound state. This is the case for sufficiently small mx since the unkown Kahler 
terms are higher order in mx, and the superpotential splitting dominates. Combining this 
with West's result suggests that in SYM the glueball states may actually be lighter than the 
S states near the SUSY point! (In other words, the S states cannot be dominating West's 
inequalities.) If so, lattice measurements of properties of the AA(a;) correlation functions will 
actually be dominated by those light glueballs (due to the mixing), and we will not obtain 
any information on the VY model. Of course, another possibility is that the glueball-S 
mixing is so strong that the VY model is not a very good description in any case. 



Finally we note, as pointed out in [|T6| , that the VY model apparently has an extra SUSY 
vacuum corresponding to (0) = 0. At this point the expectation value of S is singular and 
so it is not clear how to interpret this vacuum. Shifman and Kovner have proposed that the 
vacuum is real and represents some conformal, Z2Nc preserving point of the theory. It would 
be interesting to look for this vacuum in lattice simulations but unfortunately as can be seen 
from Eq.(pTD there is no value of soft breaking mass for which such a vacuum would be the 
global minimum. This will make it difficult to observe in lattice simulations. We return to 
a possible lattice signature of this vacuum in the following section. 



5 Tuning to SUSY 

In this section we discuss the problem of tuning a lattice simulation toward the SUSY limit. 
As in the case of chiral symmetry in QCD, the lattice regularization explicitly violates SUSY 
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[^. Gauge invariance only allows a single renormalizable operator that breaks SUSY, the 
gaugino mass. One must tune the bare gaugino mass to a special value in order to ensure 
that the symmetry (in this case SUSY) is restored, and even then it is only fully recovered in 
the continuum limit. It is important to identify a sensitive diagnostic with which to fine-tune 
the bare gaugino mass. (Of course the mass splittings described in the previous section are a 
possibility, but they are comparatively slowly varying with mx.) This is particularly true, as 
we will discuss below, because one of the most interesting quantities we wish to compute is 
the gaugino condensate (AA) which depends sensitively on small SUSY-breaking effects. We 
are tempted to propose the partition function itself as a diagnostic, since it must approach 
unity in the SUSY limit, where the vacuum energy, eo, is zero, 

Z = exp[-Veo] 1 . (25) 

The behavior of the vacuum energy is extremely sensitive to the SUSY limit. In a 
continuum calculation, one obtains an expression of the form 

eo = + m^a-^ + mA^ + A^C((Aa)'=) . (26) 

Here m is the effective gaugino mass, which vanishes when the bare gaugino mass mx is 
properly chosen. When the gaugino mass is non-zero, the vacuum energy receives quadrat- 
ically divergent corrections, resulting from non-cancellation of gaugino and gauge vacuum 
loops. Since the (bare) gaugino condensate can be related to the derivative of Z with respect 
to mx, it too is very sensitive to SUSY breaking corrections, which will appear multiplied 
by two powers of the UV cutoff. If m is not tuned sufficiently close to zero, the resulting 
contamination will preclude a measurement of the non-zero condensate which remains in the 
SUSY limit. 

Unfortunately, due to the nature of the lattice regularization, the vacuum energy eo 
does not actually approach zero, even in the SUSY limit. This is because SUSY can only 
be recovered in the IR and never in the UV part of the lattice model, which nevertheless 
contributes to eo. This is easy to see, because the lattice dispersion relations of fermions and 
bosons differ significantly in higher orders of (ka). SUSY cannot hold even approximately 
for modes with momentum k ~ a~^. 

A better method for identifying the point in parameter space corresponding to super- 
symmetry has been proposed by Montvay [0. As can be seen from Eq.(^) the shift in 
energy of the Nc SYM vacua depend on the phase of the gaugino mass in the softly broken 
theory. The lattice simulations are restricted to real mass terms in order to allow Monte 
Carlo techniques to be employed but the mass can be tuned through zero to negative values. 
For even values of Nc the n = and n = Nc/2 vacua interchange as the true minima when 
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the sign of m\ is flipped. For odd Nc the n = vacua is prefered for positive masses, while 
a negative real mass places the system on the edge of a first order phase transition between 
two vacua with conjugate phases. 

We can use the phase transition between the different Zn^ vacua as a rough indicator of 
SUSY. This transition occurs near m = (i.e. as the mass term switches sign), although 
there will be a slight overshoot due to supercooling. (For \m\ sufficiently small the critical 
bubble necessary for the transition will be larger than the lattice volume.) Note that the sign 
of the overshoot depends on the direction from which m = is approached (hysteresis)[]. We 
can define the corresponding phase transition points as mx{±), and average them to obtain 
the true SUSY point: 

ml = {mx{+) + m,(-))/2, (27) 

where 

m(m^) = 0. 

In practice, mx{±) should be defined as the points where some specified order parameter 
deviates by some specified amount from its behavior in the pure phase. A possible order 
parameter is the Wilson loop, (W), which is almost independent of m for small m, and hence 
will only display jumps near the transition points. 

It is possible that the existence of the extra vacuum proposed in and discussed above 



can be ascertained by the behavior of the Wilson loop near the phase transition. Consider 
the average of a large Wilson loop: (W). We expect that, due to confinement, this exhibits 
area law behavior with some string tension k: 

(W) kA + ■■■ , 

where the ellipsis indicate subleading effects which scale like the perimeter of the loop. 
Note that since the gluinos are in the adoint representation they, like the gluons, cannot 
fully screen sources in the fundamental representation, which must be connected by gluonic 
strings. 

Now consider what happens to (W) at the phase transition, near m = 0. If the system 

only exhibits the two vacua with non-zero gluino condensate ( (A A) = ± 327r^A^ ), then 

the important configurations near the transition will consist of regions of each of these 

phases, with domain walls in between. In the bulk of each region the string tension will be 

essentially k, assuming that the supercooling is negligible and |m| is very small (this requires 

a large lattice). On each Euclidean time slice of the loop, one can imagine the gluonic string 

^For the system to exhibit hysteresis, some "memory" of the path in m is required. This would be the 
case if, as m is varied, the previous configurations are used as the initial ones for the subsequent Monte 
Carlo update. This allows the system to be trapped in a metastable phase for mx sufficiently close to m^. 
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connecting the sources at either end, but passing through domain walls and regions of each 
phase in between. The area law part of (W) will remain the same even at the transition, 
although the subleading perimeter effects may exhibit a discontinuity due to the sudden 
appearance of the domain walls. Therefore we do NOT expect a discontinuity in the leading 
behavior of the Wilson loop. 



On the other hand, if a phase exists with (AA) = 0, as has been suggested in [0, it 
would presumably have very distinct properties including a string tension n' ^ k. In this 
case, at the phase transition the dominant configurations will contain all three phases, and 
the effective string tension will be altered to some value in between k and k'. This would 
manifest itself as a very rapid change in the leading behavior of (VT). An alternative method 
of searching for the (AA) = vacuum is to use the "multicanonical" method, which induces 
transitions between the different vacua. If the simulation is tuned sufficiently close to the 
SUSY point, it will then spend a significant portion of its time in the exotic phase. 

Also worth investigation is the difference in vacuum energy between the metastable and 
stable vauca, defined by 

[\V{my) = er(^A) - e'oimx) . (28) 

AV is finite, and vanishes at the exact SUSY point. The behavior of AV{m\) is predicted 
by (|T0|) , which we emphasize follows directly from the (Seiberg) effective Lagrangian (1), 
independent of the further V-Y analysis concerning the bound states' masses Further- 



more, the slope of AV is directly related to the gaugino condensate. Alternatively, one could 
also directly compute the quantity 

hm [ (AA)'"^(m,) - (AA)^(m,) ] , (29) 

in which the divergences also cancel. A sufficiently accurate measurement would also provide 
a determination of a, as can be seen from Eg . ([TT|) . 



Finally, we mention the possibility of measuring the surface tension of domain walls 
separating two different phases. A domain wall configuration can be produced by splitting 
the lattice into two separate regions, each treated with different (say, opposite) values of mx- 
Recently, Shifman and collaborators have given exact expressions for the profiles and energy 



densities of such configurations, using the = 1 SUSY algebra with central extension |18 



A subtlety in the preceding discusion is that our expressions are given in terms of the bare 



quantities in the "Wilsonian" regularization scheme (e.g. DR [T^), whereas in the lattice 
simulations it is the bare lattice parameters which are varied. Fortunately, one can relate 
both quantities to the bare parameters in the usual MS scheme using perturbation theory. 
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6 Summary and Prospects 



To summarize, we believe the following tests of the (softly broken) SUSY results will be 
feasible: 

(1) The first essential step is to confirm the phase structure of the theory. As discussed in 
the introduction it is possible that the SUSY phase ceases to exist for any non-zero gaugino 
mass, although one expects that it will persist until some critical value of m, possibly m ~ A 
if the gaugino plays an important role in the dynamics of the SUSY theory. Assuming that 
the SUSY phase persists, it should be described by the softly broken theory and one would 
hope to test the Zn^ vacuum structure. This structure implies the existence of a phase 
transition near the SUSY point identified above. There is even the possibilty of using the 
Wilson loop (W) as a probe of Shifman's putative vacuum at zero gaugino condensate |16[ . 



(2) Measurement of the gaugino condensate and fundamental scale A. This would provide a 
direct test of Seiberg's tower of SQCD superpotentials. The shift in the gaugino condensate 
with increasing mx could also be measured and would determine the parameter a. 

(3) Measurement of mass-splitting ratios rh (0), which are independent of the relationship 
between different regularizations. This is a test of the V-Y Lagrangian itself. At the SUSY 
point, the masses should be degenerate and provide a measurement of the product aA. 

(4) On a more speculative note, some of the recent results of Shifman and collaborators 
[T^ are potentially testable, including the controversial existence of a Z2Nc^pi'eserving 

vacuum and the exact calculation of domain wall energy densities. 

Finally, we mention a related theory which has been speculated to have very different 
behavior from normal QCD and which could be easily simulated using current lattice tech- 
nology. The model is an SU(2) gauge theory with a Dirac spinor of adjoint fermions. This 
model was studied in |jT9| with real Susskind fermions in the strong coupling expansion where 



it gives rise to a massless composite Dirac fermion. This theory can be reached at tree level 
by the inclusion of a large soft breaking scalar mass in N=2 SQCD, although in that case 



the spur ion symmetries are not sufficient to determine the IR behavior, so the results of 119 
cannot be confirmed without further analysis. 
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7 Appendix: Rescaling Anomaly 



The results given in this paper imphcitly assumed the holomorphic normahzations of fields 
in the bare Lagrangian. A rescaling is necessary in order to compare our predictions to 
lattice results which result from canonically normalized bare Lagrangians (i.e. unit kinetic 
energy coefficients for the gluino, as opposed to 1/(7^). Due to the rescaling anomaly, some 
care is necessary in this rescaling . 



The vector superfields in the holomorphic and canonical normalizations are related by 

Vh = gcVc , (1) 

which, in particular, implies Xh = gc^c- Naively, using classical rescaling, we would have the 
following relation between the condensates: 

(A,A,) = gliKK) . (2) 

However, the rescaling anomaly introduces additional effects. In the holomorphic computa- 
tion, we essentially computed 



J DVh XhXh exp (-J^/ d^^d^^ WhWh 



(3) 



Changing variables to Vc, we have 

J Dig^V,) gXX, exp (^-(^ -^\ng,)^J d'xdH WM^j , (4) 



where the shift in coupling by ^ In gc is due to the rescaling anomaly, arising from the 
functional measure. We see that if we take 

the path integral in is canonically normalized, with coupling constant gc- 
Our prediction for the canonically normalized lattice calculation is therefore 

(AA) lattice.gc = -^3271^ A3, (6) 

gl 

where A is the holomorphic strong couping scale. A is defined by the RGE evolution of gh, 
governed by the exact one loop beta function, and is distinct from the strong coupling scale 
associated with gc- 
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